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Abstract 

For single- j shells with j = |, § and we relate the large probability of 
I + ground states to the largest (smallest) coefficients Qj/^m = (nvf3I\ A J ^ • 
A J \nv(3I), where n is the particle number, v is the seniority, (5 is an additional 
quantum number, and I is the angular momentum of the state. Interesting 
regularities of the probabilities of I + ground states are noticed and discussed 
for 4-particle systems. Several counter examples of the + ground state (OGS) 
predominance are noticed for the first time. 
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The shell structure and pairing are "universals" which play important roles in many- 
body systems such as metal clusters, atoms, and nuclei etc. Therefore, a study which relates 
intrinsic properties of the shell model with the origin of the pairing phenomenon in many- 
body systems is extremely important. Recent work by Johnson et al |T]-§] was the first 
effort along these lines. In that work, the low-lying spectra of many-body systems with 
even number of particles were examined using random two-body interactions (TBRE), and 
its results showed a preponderance of I w = + ground states (OGS). This phenomenon also 
occurs in s<i-boson systems, and was referred to as "a robust property" of the valence space 
All these studies demonstrate that the features of pairing arise from a very large 
ensemble of two-body interactions and seem to be independent of the specific character of 
the force. These observations, however, seem to be contrary to the traditional assumption 
in nuclear physics, whereby the OGS dominance in even-even nuclei is a reflection of a 
strong pairing interaction associated with a strong short-range attraction between identical 
nucleons. 

It is therefore very important and interesting to study why the OGS are favored in even- 
even nuclei when using the TBRE There have been several efforts to understand 
this phenomenon. In Ref. |J an analysis was carried out on systems of identical particles 
occupying orbits with j = |, |, and |, and it was indicated that there is a correspondence 
between a larger distribution width of + states and the OGS dominance. In another paper, 
it was proposed ]7| that for a system of interacting bosons the probability that the ground 
state has a certain value of the angular momentum is not really fixed by the full distribution 
of eigenvalues, but rather by that of the lowest one. In ||, Muhall et al. discussed the 
OGS dominance within single- j shells by using geometric chaoticity and uniformly changed 
random interactions. In ||, Kusnesov discussed the sp-boson case. In [jn| Kaplan et al. 
studied the correlation between the eigenvalues and spins of the states in some simple cases. 

All these studies [1-10] on the OGS dominance are interesting and important, and have 
potentially impacted our understanding on the origin of one of the most characteristic fea- 
tures of nuclear spectra. An essential understanding of the OGS predominance, however, has 
not yet been achieved. 

Towards that goal, we present in this Rapid Communication a new view on the origin of 
the OGS dominance. We are able to achieve this new understanding because we focus on the 
simplest non-trivial problem, that of identical nucleons in a single- j shell. This problem has 
the advantage that it is simply solvable, thus permitting useful insight to be obtained which 
can hopefully then be applied to more general problems. In the process, some previously 
unrecognized features, such as counter examples of the OGS predominance, and regularities 
of probability of I + ground states, have been discerned. 

We first discuss the properties of linear combinations of random numbers. Let G j be a 
set of random numbers with a distribution function 

p(Gj) = -^exp(-G 2 j), J = 0, 2, • • ■ 2(j - 1), (0.1) 
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and F(k) be a set of linear combinations of G 
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F{k) = Y,<*iGj, k = l,2,---N, (0.2) 
j 

where j labels the single- j shell, and N is the total number of states (to be introduced later). 
It can be shown that distribution functions of random F(k) are 

p{F{k)) = -^=exp(- (F(k)f /g 2 k ), g\ = £ (a{f . (0.3) 
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If in Eq. ( p.2|) is the largest (or the smallest) among all the a{ (k = l,---N), the 
probability of F(m) being both the smallest and the largest number is large. To show this, 
let us look at 

F(k) = F{k) - F{m) = c J k 'Gj, + ( £ (a{ - a J m ) Gj) = c J k 'Gj, + F(k), (0.4) 

where k ^ m, and c-f = (a k — a^), J~'(k) = — «m) @J- The right hand side of 

Eq. (p.4[) has two terms, both of which are Gaussian type random numbers. The coefficient 
qf is negative or positive, and thus effectively produces a shift in T{k), as is evident from 
Eq. Q0.4p . Therefore, all of the functions J-{k) have large probabilities to be both negative 
and positive, depending on the sign of the shift, i.e., F(m) has a large probability to be both 
the smallest and the largest. 

Roughly speaking, one may use d — \c{',\, which is the smallest among all the \cj. \(k ^ m), 

as an estimate of the shifts, and use the maximum of D = \Jy^j^j' ( a k ~ a mf as an estimate 
of the distribution widths of the functions T'{k). The probability of F(m) being the smallest 
(largest) is very large when d is comparably large (d ~ D). This occurs when there are a 
very few (3~6) parameters in the Hamiltonian. The probability of F(m) being the smallest 
(largest) is roughly determined by g k , the widths of the quantities F(k), when d << D. The 
states with large eigenvalue widths have large probabilities to be both the smallest and the 
largest. This situation occurs if there are many (e.g., more than 20) TBRE parameters. In 
Ref. 0, the OGS probability are determined by the total widths of the energy eigenvalues 
because there are 30 parametrized T=l interactions. 

If there are two or more coefficients (J' = 0, 2, • • - 2j — 1) which are the largest or 
smallest for different functions F(k), the probability of finding F(m) as both the smallest 
and the largest increases. 

Keeping the above ideas in mind, we now return to the shell model calculation for 
single— j shells. The Hamiltonian for a single— j shell is 



H = V2J+lGj [A J 1 x i J 
j 

A J t = -L [ a J x a^ J ,A J = --L [a 3 x a/ ,Gj = (j 2 J\V\j 2 J). 
The parameters Gj are distributed according to Eq. ( p.l| ). 
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Consider the simplest example, a j = 7/2 shell with 4 particles, where all the states are 
labeled by their total angular momenta / and their seniority quantum numbers (v). The 
eigenvalues are as follows PTTf] : 

Eo(0) — + -G2 + -G4 + ~~^~Gq, E 2 (2) = -Go + ~^G 2 + -G4 + ~^Gq, 

-^2(4) = G 2 + + YJ^ 6 ' -^ 4 ( 2 ) = g^ 2 1T^ 6 ' 

-^4(4) = -C 2 + 1 ■ G 4 + -G 6 , £5(4) = -G 2 + ~^-G± + —G 6 , 

E 6 (2) = -G + ~G 2 + -G4 + —G 6 , E 8(4 ) = ^-^2 + ~^f G ^ + ~33" G6 ' (°' 5 ) 

where bold font is used for the largest and italic for the smallest amplitudes in an expansion 
in terms of Gj, 

E I{v) =^ J m G,j. (0.6) 
j 

Using Eq. (p.5|), we can predict the probability of I + ground states without carrying out 
energy calculations with the TBRE. For example, the OGS probability is determined by 



(IGq / (IG2 / dGi I cLGq I dEo(o) / dE 2 ( 2 ) ' ' ' / ^8(4) 

J J J J ■'#0(0) ^#0(0) 

* (^o(o) - £ «o J (0 )G^ • • ■ 5 ^ 8 (4) - £ a 8 J (4) G^ p(G )p(G 2 )p(G 4 )p(G 6 ). (0.7) 

The probabilities of finding / as the spin of ground states for a 4-particle system with j — I 
are listed in Table I. The row "test" corresponds to results obtained by using Eq. (P.5|) and 
1000 sets of the TBRE. The row "cal." corresponds to the probabilities predicted by an 
integral for each I + state similar to Eq. ( p.7|) for the + state. The probabilities calculated 
using the TBRE and those predicted by using integrals like Eq. (p.7|) are all consistent with 
one another. 

Using Eq. (|0.5|) and Eq. (p.3|), gi( v )— the distribution width of Ei( v ) — is obtained and 
listed in the last row of Table I. It is found that there is no correspondence between the 
width gi( v ) and the probability of finding I + (v) to be the ground state. 

Because there are only 4 parameters in Eq. (|0.5| ) , we can use the shift argument to predict 
which states have large probabilities to be the ground. They are the states with the largest 
(smallest) ocj {v) : l(v)=0(0), 2(4), 4(4), 8(4). 

Our success in explaining which states can have large probabilities to be the ground state 
for the j — \ shell encouraged us to go to higher j shells. The study of 4 particles in the 
j = I shell is very interesting because the OGS is predominant in this case (See Fig. 1). 
This is therefore a good case in which to check whether the above assumptions continue to 
apply. 



4 



Unfortunately, the eigenstates of 4 particles in this shell are complicated by mixings 
between states with the same and different seniorities. In such a case there is no simple 
relation between the eigenvalues and the two-body interaction parameters, Gj. 

To simplify the analysis for the j = § case, we assume that the mixings are negligible, and 



(nv 



pi\EjV2J + T 
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j\\J x j^J \nvf3I) are the "eigenvalues". This assumption is supported 
by the fact that in the basis \nvf3I) |T2[ mixings are observed to be small for / > states 
in even-iV systems in the j = I shell. It is also supported by a recent claim [§] that a 
many-body system picks out seniority as an approximate quantum number, even though it 
is not obviously implicit in the Hamiltonian. If we could take those mixings into account 
in a simple way, those off-diagonal matrix elements would presumably pull down the lower 
"eigenstates" even lower, and push up the higher ones even higher. For 1 = states, the 
off-diagonal matrix elements are comparable in magnitude to the diagonal matrix elements, 
and the pulling-down and pushing-up effects should be more important than for 1^0 states. 
Therefore, this assumption leads to an underestimate of the OGS probability. 

The coefficients aj for each Gj in the "eigenvalues" are listed in Table II. There are 
three largest aj (J = 0, 4, 6) and one smallest aj (J = 8)in the 1 = states. Therefore, the 
OGS probability must be very large. In addition, there are one largest aj (J = 2) and one 
smallest aj (J = 6) in one of the 4 + states. Finally, there are one largest aj (J = 8) and 
two smallest aj (J = 2,4) in the highest angular momentum state (J = 12). This indicates 
large probabilities of 4 + ground states and 12 + ground states. In one of the 2 + states (v = 2) 
the coefficient of Gq is very near to the smallest. Thus, there should be a sizable percentage 
of 2 + ground states as well. 

Shell model calculations based on the TBRE are consistent with the above discussion: 
The probability of + states to be the ground state is 66.4%, that of 4 + is 12.2% and that 
of 12 + is 17.4%. The probability of finding 2 + states as the ground state is 3.4%. There 
are very small probabilities (0.6% in total) for all other I + ground states. The predicted 
OGS probability obtained using an integral similar to Eq. (p.7|) , where no mixing is assumed 
(which is exact for j = |), is 45.11%. As expected, this is smaller than the probability 
(66.4%) found when the mixings are taken into account. We note without further details 
that the above argument is also applicable to 4 and 6 particles in the j = y shell, and to 
the general features of odd-A fermion system in small single- j shells and srf-boson systems 

We have also carried out several other sets of calculations, both for larger single- j shells 
and in some cases for two-j shells. Some of the principal outcomes are: 

1) The probability of ground states with odd angular momenta is much smaller than 
that of their neighboring even angular momenta even if the number of states with the same 
angular momentum I found in j n configuration is comparably large. 

2) The unique I max state has a large probability to be the ground state, although this 
probability decreases with j. According to ref. ||, by using random interactions which 
distribute uniformly between —1 and 1, the probability of I^ ax ground states staggers rapidly 



5 



and becomes for several single j. 

3) The I = 2 + and 4 + states have large probabilities to be the ground state. This 
indicates that small and even angular momentum states are favored as the ground state in 
single- j shells. 

4) The OGS predominance obtained using the TBRE is not a "rule" without exceptions. 
There are counter examples: j ' = | and j = y sne U s with 4 particles (refer to Fig. 1); two-j 
shell (j = | plus |) with 4 particles, where the probability of the OGS is 24% while that of 
the 2+ GS is 33%. 

A summary of the above regularities is shown in Fig. 1, from which we see that when 
the OGS probability is a maximum, the probability of 2 + ground states is a minimum, and 
vice versa. Fig. 1 indicates that the fluctuation of the OGS probability decreases with j, the 
OGS probability saturates around 35%, and the GS probabilities of / = 2,4 and I max states 
decrease and converge around 9%. These regularities are worthy of further study. 

To summarize, we have presented in this Rapid Communication an interpretation of the 
OGS probability in terms of the largest and smallest coefficients cnjr v -\ for several single- j 
shells. The interpretation applies when seniority is either a conserved or an approximately 
conserved quantum number. Further study showed that this method is not only applicable 
to even number of fermions in a small single- j shell, but also is a good benchmark to explain 
the /GS probability of odd-A fermions in a small single- j shell and general features of sd- 
boson systems |12|]. Therefore, for the first time, the JGS probabilities of a variety of systems 
can be discussed on the same footing using the method proposed in this paper. 

The OGS probabilities have been calculated using integrals like Eq. (|0.7| ), the predicted 
probabilities being reasonably consistent with those obtained using the TBRE. We showed 
that the width of the distribution of + eigenvalues is not a manifestation of the OGS 
preponderance in small single- j shells. 

Two other conclusions not yet mentioned are as follows: 1) If I + states have a very large 
probability to be the ground state, there must be, on average, a large energy gap between 
the I + ground states and other states. This gap comes from a large shift d in small j shells, 
where there are very few parameters, or from the difference between the width of eigenvalues 
for the I + states and those for all other angular momentum states in very large single- j shells 
and multi-j shells. 2) If certain I + states are favored to be the ground state, they are also 
favored to be the highest state with the same probability. 

As j increases, the approximations made in this study, i.e., the omission of mixings 
between states with different seniorities and those among states with the same seniority 
deteriorate. In principle, the relationship between the eigenvalues and the two-body inter- 
actions is not linear. As a consequence, more general or more complicated cases may not 
be studied in this simple way. Nevertheless, our analysis should provide helpful clues for 
what is going on with regards to OGS predominance in more general even-even many-body 
systems. 

Finally, it is not yet understood why certain coefficients turn out to be the largest 



6 



and/or the smallest among all of them. Since this observation seems to be critical to ex- 
plaining OGS dominance, at least in small single- j shells, further consideration of this issue 
is warranted. 

We are grateful to Dr. N.Yoshinaga for his generosity in allowing us to use his modified 
"two-j shell" program. Many discussions with Drs. N. Yoshinaga, S. Pittel, and S. Yamaji, 
and communications with Dr. C.W. Johnson are gratefully acknowledged. We thank Dr. 
S. Pittel for his careful reading of the mansucript. This work is supported in part by the 
Japan Society for the Promotion of Science (contract ID: P0102) and the U.S. NSF under 
grant No. PHY-9970749. 
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Caption: 

Fig. 1. Probabilities of I + ground states for different j shells with 4 particles. All 
probabilities are obtained from 1000 runs of the TBRE. 
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TABLE I. The probability of each state to be the ground state and its distribution width 
for 4 particles in a j = 7/2 shell. All eigenstates are uniquely labeled by their seniorities 
and angular momenta. The first line specifies the angular momentum and seniority for each 
state. The probabilities labeled by "test" are obtained from 1000 runs of the TBRE, and 
those labeled by "cal." are obtained by integrals like Eq. ( p.7| ) for + state (see the text for 
further details). The distribution width, <?/(„), of each energy level is listed in the 3rd row. 
The shift, d, of each level is given in the last row. The d inside a bracket corresponds to a 
shift coming from the smallest ctj/ v y 



I(v) 


0(0) 2(2) 


2(4) 


4(2) 


4(4) 


5(4) 6(2) 


8(4) 


test 


19.9% 1.2% 


31.7% 


0.0% 


25.0% 


0.0% 0.0% 


22.2% 


cal. 


18.19% 0.89% 


33.25% 


0.00% 


22.96% 


0.00% 0.02% 


24.15% 


9i(v) 


3.14 3.25 


4.12 


3.45 


3.68 


3.62 3.64 


4.22 


d 


1.00 


1.32 (0.98; 


) - 


0.50(0.50) 




0.68(0.35) 



TABLE II. The coefficients &i( v p\ for 4 particles in a j — | shell. Bold font is used for 
the largest aj are the largest and italic for the smallest aj. 



I 


v 


Go 


G 2 


G4 


Gq 


Gg 








1.60 


0.50 


0.90 


1.30 


1.70 





4 





.00 


0.20 


2.57 


2.91 


0.32 


2 


2 





.60 


1.43 


1.22 


0.893 


1.86 


2 


4 





.00 


1.35 


1.69 


1.70 


1.26 


3 


4 





.00 


0.36 


2.28 


2.63 


0.71 


4 


4 





.00 


0.50 


2.08 


2.43 


0.99 


4 


4 





.00 


2.04 


1.02 


0.890 


2.06 


4 


2 





.60 


0.68 


1.04 


2.40 


1.28 


5 


4 





.00 


1.00 


1.59 


1.84 


1.57 


6 


4 





.00 


1.64 


0.98 


1.08 


2.29 


6 


4 





.00 


0.39 


1.85 


2.34 


1.43 


6 


2 





.60 


0.34 


1.66 


1.33 


2.07 


7 


4 





.00 


1.20 


1.09 


1.40 


2.31 


8 


4 





.00 


0.41 


1.42 


2.05 


2.13 
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,60 


0.55 


0.68 


1.58 


2.59 
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4 





,00 


0.17 


1.33 


2.12 


2.38 


10 4 





.00 


0.70 


0.69 


1.41 


3.21 


12 4 





.00 


0.00 


0.52 


1.69 


3.78 
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